Molecular Evidence for Relaxed Selection on the Enamel Genes of Toothed Whales (Odontoceti) with Degenerative Enamel Phenotypes.

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development.

Drivers of morphological evolution in the toothed whale jaw.

Toothed whales (odontocetes) emit high-frequency underwater sounds (echolocate)-an extreme and unique innovation allowing them to sense their prey and environment. Their highly specialized mandible (lower jaw) allows high-frequency sounds to be transmitted back to the inner ear. Echolocation is evident in the earliest toothed whales, but little research has focused on the evolution of mandibular form regarding this unique adaptation. Here, we use a high-density, three-dimensional geometric morphometric analysis of 100 living and extinct cetacean species spanning their ∼50-million-year evolutionary history. Our analyses demonstrate that most shape variation is found in the relative length of the jaw and the mandibular symphysis. The greatest morphological diversity was obtained during two periods of rapid evolution: the initial evolution of archaeocetes (stem whales) in the early to mid-Eocene as they adapted to an aquatic lifestyle, representing one of the most extreme adaptive transitions known, and later on in the mid-Oligocene odontocetes as they became increasingly specialized for a range of diets facilitated by increasingly refined echolocation. Low disparity in the posterior mandible suggests the shape of the acoustic window, which receives sound, has remained conservative since the advent of directional hearing in the aquatic archaeocetes, even as the earliest odontocetes began to receive sounds from echolocation. Diet, echolocation, feeding method, and dentition type strongly influence mandible shape. Unlike in the toothed whale cranium, we found no significant asymmetry in the mandible. We suggest that a combination of refined echolocation and associated dietary specializations have driven morphology and disparity in the toothed whale mandible.

Speciation in the deep: genomics and morphology reveal a new species of beaked whale Mesoplodon eueu.

The deep sea has been described as the last major ecological frontier, as much of its biodiversity is yet to be discovered and described. Beaked whales (ziphiids) are among the most visible inhabitants of the deep sea, due to their large size and worldwide distribution, and their taxonomic diversity and much about their natural history remain poorly understood. We combine genomic and morphometric analyses to reveal a new Southern Hemisphere ziphiid species, Ramari's beaked whale, Mesoplodon eueu, whose name is linked to the Indigenous peoples of the lands from which the species holotype and paratypes were recovered. Mitogenome and ddRAD-derived phylogenies demonstrate reciprocally monophyletic divergence between M. eueu and True's beaked whale (M. mirus) from the North Atlantic, with which it was previously subsumed. Morphometric analyses of skulls also distinguish the two species. A time-calibrated mitogenome phylogeny and analysis of two nuclear genomes indicate divergence began circa 2 million years ago (Ma), with geneflow ceasing 0.35-0.55 Ma. This is an example of how deep sea biodiversity can be unravelled through increasing international collaboration and genome sequencing of archival specimens. Our consultation and involvement with Indigenous peoples offers a model for broadening the cultural scope of the scientific naming process.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos.

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae).

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species.

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

Positive Selection and Inactivation in the Vision and Hearing Genes of Cetaceans.

The transition to an aquatic lifestyle in cetaceans (whales and dolphins) resulted in a radical transformation in their sensory systems. Toothed whales acquired specialized high-frequency hearing tied to the evolution of echolocation, whereas baleen whales evolved low-frequency hearing. More generally, all cetaceans show adaptations for hearing and seeing underwater. To determine the extent to which these phenotypic changes have been driven by molecular adaptation, we performed large-scale targeted sequence capture of 179 sensory genes across the Cetacea, incorporating up to 54 cetacean species from all major clades as well as their closest relatives, the hippopotamuses. We screened for positive selection in 167 loci related to vision and hearing and found that the diversification of cetaceans has been accompanied by pervasive molecular adaptations in both sets of genes, including several loci implicated in nonsyndromic hearing loss. Despite these findings, however, we found no direct evidence of positive selection at the base of odontocetes coinciding with the origin of echolocation, as found in studies examining fewer taxa. By using contingency tables incorporating taxon- and gene-based controls, we show that, although numbers of positively selected hearing and nonsyndromic hearing loss genes are disproportionately high in cetaceans, counts of vision genes do not differ significantly from expected values. Alongside these adaptive changes, we find increased evidence of pseudogenization of genes involved in cone-mediated vision in mysticetes and deep-diving odontocetes.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture.

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for $\sim $3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from $>$38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.].

The complete mitochondrial genome of the critically endangered Atlantic humpback dolphin, Sousa teuszii (Kükenthal, 1892).

The Atlantic humpback dolphin remains an understudied, critically endangered cetacean species. Here, we describe the first complete mitogenome of Sousa teuszii, derived from an animal stranded on Île des Oiseaux, Sine Saloum, Senegal. The S. teuszii mitogenome is composed of 16,384 base pairs and is 98.1% identical to its closest relative with a mitogenome, Sousa chinensis. Phylogenetic analysis confirms its placement with S. chinensis, as well as the placement of the genus Sousa within subfamily Delphininae.

The core transcriptome of mammalian placentas and the divergence of expression with placental shape.

INTRODUCTION: The placenta is arguably the most anatomically variable organ in mammals even though its primary function is conserved. METHOD: Using RNA-Seq, we measured the expression profiles of 55 term placentas of 14 species of mammals representing all major eutherian superordinal clades and marsupials, and compared the evolution of expression across clades. RESULTS: We identified a set of 115 core genes which is expressed (FPKM ≥10) in all eutherian placentas, including genes with immune-modulating properties (ANXA2, ANXA1, S100A11, S100A10, and LGALS1), cell-cell interactions (LAMC1, LUM, and LGALS1), invasion (GRB2 and RALB) and syncytialization (ANXA5 and ANXA1). We also identified multiple pre-eclampsia associated genes which are differentially expressed in Homo sapiens when compared to the other 13 species. Multiple genes are significantly associated with placenta morphology, including EREG and WNT5A which are both associated with placental shape. DISCUSSION: 115 genes are important for the core functions of the placenta in all eutherian species analyzed. The molecular functions and pathways enriched in the core placenta align with the evolutionarily conserved functionality of the placenta.

Ancestral resurrection of anthropoid estrogen receptor ß demonstrates functional consequences of positive selection.

Anthropoid primates arose during the Eocene approximately 55 million years ago (mya), and extant anthropoids share a most recent common ancestor ∼40mya. Paleontology has been very successful at describing the morphological phenotypes of extinct anthropoids. Less well understood is the molecular biology of these extinct species as well as the phenotypic consequences of evolutionary variation in their genomes. Here we resurrect the most recent common ancestral anthropoid estrogen receptor ß gene (ESR2) and demonstrate that the function of this ancestral estrogen receptor has been maintained during human descent but was altered during early New World monkey (NWM) evolution by becoming a more potent transcriptional activator. We tested hypotheses of adaptive evolution in the protein coding sequences of ESR2, and determined that ESR2 evolved via episodic positive selection on the NWM stem lineage. We separately co-transfected ESR2 constructs for human, NWM, and the anthropoid ancestor along with reporter gene vectors and performed hormone binding dose response experiments that measure transactivation activity. We found the transactivation potentials of the ancestral and human sequences to be significantly lower (p<0.0001 in each comparison) than that of the NWM when treated with estradiol, the most prevalent estrogen. We conclude the difference in fold activation is due to positive selection in the NWM ERß ligand binding domain. Our study validates inferential methods for detecting adaptive evolution that predict functional consequences of nucleotide substitutions and points a way toward examining the functional consequences of positive Darwinian selection.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals.

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

A phylogenomic analysis of the role and timing of molecular adaptation in the aquatic transition of cetartiodactyl mammals.

Recent studies have reported multiple cases of molecular adaptation in cetaceans related to their aquatic abilities. However, none of these has included the hippopotamus, precluding an understanding of whether molecular adaptations in cetaceans occurred before or after they split from their semi-aquatic sister taxa. Here, we obtained new transcriptomes from the hippopotamus and humpback whale, and analysed these together with available data from eight other cetaceans. We identified more than 11 000 orthologous genes and compiled a genome-wide dataset of 6845 coding DNA sequences among 23 mammals, to our knowledge the largest phylogenomic dataset to date for cetaceans. We found positive selection in nine genes on the branch leading to the common ancestor of hippopotamus and whales, and 461 genes in cetaceans compared to 64 in hippopotamus. Functional annotation revealed adaptations in diverse processes, including lipid metabolism, hypoxia, muscle and brain function. By combining these findings with data on protein-protein interactions, we found evidence suggesting clustering among gene products relating to nervous and muscular systems in cetaceans. We found little support for shared ancestral adaptations in the two taxa; most molecular adaptations in extant cetaceans occurred after their split with hippopotamids.

The evolution of embryo implantation.

Embryo implantation varies widely in placental mammals. We review this variation in mammals with a special focus on two features: the depth of implantation and embryonic diapause. We discuss the two major types of implantation depth, superficial and interstitial, and map this character on a well-resolved molecular phylogenetic tree of placental mammals. We infer that relatively deep interstitial implantation has independently evolved at least eight times within placental mammals. Moreover, the superficial type of implantation represents the ancestral state for placental mammals. In addition, we review the genes involved in various phases of implantation, and suggest a future direction in investigating the molecular evolution of implantation-related genes.

Genome-wide adaptive complexes to underground stresses in blind mole rats Spalax.

The blind mole rat (BMR), Spalax galili, is an excellent model for studying mammalian adaptation to life underground and medical applications. The BMR spends its entire life underground, protecting itself from predators and climatic fluctuations while challenging it with multiple stressors such as darkness, hypoxia, hypercapnia, energetics and high pathonecity. Here we sequence and analyse the BMR genome and transcriptome, highlighting the possible genomic adaptive responses to the underground stressors. Our results show high rates of RNA/DNA editing, reduced chromosome rearrangements, an over-representation of short interspersed elements (SINEs) probably linked to hypoxia tolerance, degeneration of vision and progression of photoperiodic perception, tolerance to hypercapnia and hypoxia and resistance to cancer. The remarkable traits of the BMR, together with its genomic and transcriptomic information, enhance our understanding of adaptation to extreme environments and will enable the utilization of BMR models for biomedical research in the fight against cancer, stroke and cardiovascular diseases.

Molecular evolution tracks macroevolutionary transitions in Cetacea.

Cetacea (whales, dolphins, and porpoises) is a model group for investigating the molecular signature of macroevolutionary transitions. Recent research has begun to reveal the molecular underpinnings of the remarkable anatomical and behavioral transformation in this clade. This shift from terrestrial to aquatic environments is arguably the best-understood major morphological transition in vertebrate evolution. The ancestral body plan and physiology were extensively modified and, in many cases, these crucial changes are recorded in cetacean genomes. Recent studies have highlighted cetaceans as central to understanding adaptive molecular convergence and pseudogene formation. Here, we review current research in cetacean molecular evolution and the potential of Cetacea as a model for the study of other macroevolutionary transitions from a genomic perspective.

The evolutionary history of cetacean brain and body size.

Cetaceans rival primates in brain size relative to body size and include species with the largest brains and biggest bodies to have ever evolved. Cetaceans are remarkably diverse, varying in both phenotypes by several orders of magnitude, with notable differences between the two extant suborders, Mysticeti and Odontoceti. We analyzed the evolutionary history of brain and body mass, and relative brain size measured by the encephalization quotient (EQ), using a data set of extinct and extant taxa to capture temporal variation in the mode and direction of evolution. Our results suggest that cetacean brain and body mass evolved under strong directional trends to increase through time, but decreases in EQ were widespread. Mysticetes have significantly lower EQs than odontocetes due to a shift in brain:body allometry following the divergence of the suborders, caused by rapid increases in body mass in Mysticeti and a period of body mass reduction in Odontoceti. The pattern in Cetacea contrasts with that in primates, which experienced strong trends to increase brain mass and relative brain size, but not body mass. We discuss what these analyses reveal about the convergent evolution of large brains, and highlight that until recently the most encephalized mammals were odontocetes, not primates.

Characterization of human cortical gene expression in relation to glucose utilization.

OBJECTIVES: Human brain development follows a unique pattern characterized by a prolonged period of postnatal growth and reorganization, and a postnatal peak in glucose utilization. The molecular processes underlying these developmental changes are poorly characterized. The objectives of this study were to determine developmental trajectories of gene expression and to examine the evolutionary history of genes differentially expressed as a function of age. METHODS: We used microarrays to determine age-related patterns of mRNA expression in human cerebral cortical samples ranging from infancy to adulthood. In contrast to previous developmental gene expression studies of human neocortex that relied on postmortem tissue, we measured mRNA expression from the nondiseased margins of surgically resected tissue. We used regression models designed to identify transcripts that followed significant linear or curvilinear functions of age and used population genetics techniques to examine the evolution of these genes. RESULTS: We identified 40 transcripts with significant age-related trajectories in expression. Ten genes have documented roles in nervous system development and energy metabolism, others are novel candidates in brain development. Sixteen transcripts showed similar patterns of expression, characterized by decreasing expression during childhood. Comparative genomic analyses revealed that the regulatory regions of three genes have evidence of adaptive evolution in recent human evolution. CONCLUSIONS: These findings provide evidence that a subset of genes expressed in the human cerebral cortex broadly mirror developmental patterns of cortical glucose consumption. Whether there is a causal relationship between gene expression and glucose utilization remains to be determined.

From PPROM to caul: The evolution of membrane rupture in mammals.

Rupture of the extraembryonic membranes that form the gestational sac in humans is a typical feature of human parturition. However, preterm premature rupture of membranes (PPROM) occurs in approximately 1% of pregnancies, and is a leading cause of preterm birth. Conversely, retention of an intact gestational sac during parturition in the form of a caul is a rare occurrence. Understanding the molecular and evolutionary underpinnings of these disparate phenotypes can provide insight into both normal pregnancy and PPROM. Using phylogenetic techniques we reconstructed the evolution of the gestational sac phenotype at parturition in 55 mammal species representing all major viviparous mammal groups. We infer the ancestral state in therians, eutherians, and primates, as in humans, is a ruptured gestational sac at parturition. We present evidence that intact membranes at parturition have evolved convergently in diverse mammals including horses, elephants, and bats. In order to gain insight into the molecular underpinnings of the evolution of enhanced membrane integrity we also used comparative genomics techniques to reconstruct the evolution of a subset of genes implicated in PPROM, and find that four genes (ADAMTS2, COL1A1, COL5A1, LEPRE1) show significant evidence of increased nonsynonymous rates of substitution on lineages with intact membranes as compared to those with ruptured membranes. Among these genes, we also discovered that 17 human SNPs are associated with or near amino acid replacement sites in those mammals with intact membranes. These SNPs are candidate functional variants within humans, which may play roles in both PPROM and/or the retention of the gestational sac at birth.